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ABSTRACT 

In this paper a partially reinforced cylinder containing an axial through 
crack is considered. The reinforcement is assumed to be fully bonded to the 
main cylinder. The composite cylinder is thus modelled by a nonhomogeneous 
shell having a step change in the elastic properties at the z=0 plane, z being 
the axial coordinate. Using a Reissner type transverse shear theory the 
problem is reduced to a pair of singular integral equations. In the special 
case of a crack tip touching the bimaterial interface it is shown that the 
dominant parts of the kernels of the integral equations associated with both 
membrane loading and bending of the shell reduce to the generalized Cauchy 
kernel obtained for the corresponding plane stress case. The integral equa- 
tions are solved and the stress intensity factors are given for various crack 
and shell dimensions. A bonded fiberglass reinforcement which may serve as 
a crack arrestor is used as an example. 

1 . Introduction 

In studying the failure of structures for the purpose of calculating 
the fracture mechanics parameters, a very large variety of structural compo- 
nents may locally be modelled as relatively thin-walled plates or shells. 
Nearly all "pressure boundaries" and piping as well as some important parts 
of aerospace and hydrospace structures may be cited as examples of such com- 
ponents. From a viewpoint of structural integrity two of the important 
questions one may be concerned with in this respect are the life estimate 
based on the subcritical growth of an existing flaw and the residual strength 
or the load carrying capacity of the structure based on the criticality of 
a dominant flaw. In most cases the two questions may be adequately dealt 
with by idealizing the component with a plate or a shell and the flaw with a 
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part-through or a through crack, by calculating the appropriate fracture 
mechanics parameter (e.g., the stress intensity factor), and by applying a 
suitable criterion along with the baseline characterization of the material. 

In relatively thin-walled structures the inherently three-dimensional 
crack problem is approximated by a "plate" or a "shell" problem, that is, by 
suppressing the thickness coordinate through the use of a plate or a shell 
theory. In the earlier studies of the subject the classical plate and shell 
theories were used to solve the problem (see, for example, [1] for review). 
However, particularly in problems requiring the calculation of the stress 
intensity factors, the necessity of using a higher order theory has now 
been well-established. For example, it has been shown that by using a Reissner 
type transverse shear theory [2] (and hence by satisfying the boundary condi- 
tions on the crack surfaces for all stress and moment resultants separately) 
one could obtain an asymptotic stress state around the crack tips which is 
identical to that given by the in-plane and anti -plane elasticity solutions 
(see, for example, [3], [4] and [5]). Furthermore, it has also been shown 
that in the limiting case of small crack lengths these plate and shell results 
approach that of plane elasticity not only for the internal but also for the 
edge cracks [6], [7], [8]. Other results obtained by using a transverse 
shear theory for various crack-shell geometries and loading conditions may 
be found in [9]-[12]. 

With the exception of [7], [8] and [13], in all crack studies in shells 
that appeared in literature the shell is assumed to be "infinite" in the 
sense that the interaction of the perturbation field of the crack with the 
boundaries of or with other geometric discontinuities in the shell are assumed 
to be negligible. The interaction of the stress field around the crack with 
a stress-free boundary and with a fully clamped boundary in a cylindrical 
shell was considered in [7] and [8], respectively, where the special case 
of the crack intersecting the boundary has also been studied. In [13] the 
interaction of a crack field with a circumferential line stiffener in a cylin- 
drical shell was studied by using the classical shell theory. In this paper 
we consider the somewhat more general problem of a nonhomogeneous cylindrical 
shell containing an axial crack. The problem studied is that of two relatively 
long cylinders having the same radii and different mechanical properties 
that are joined along their boundaries at a plane perpendicular to their 
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common axis (Fig. lb). The model may be used to simulate composite cylinders, 
cylinders with reinforcing layers, and homogeneous cylinders having a step 
change in thickness (Fig. la and b). Such solutions are needed or may be 
very useful in, for example, crack arrest studies in pipes and containers 
with reinforcements. 

2. General Formulation of the Problem 

The analytical problem under consideration is described in Fig. lb. 

The actual problem may arise, for example, from a reinforced shell shown in 
Fig. la. Since "shell theory" is used in formulating the problem, whatever 
the actual configuration and composition of the medium for X 2 > 0 » it has to 
be reduced to a homogeneous shell having the same radii as the semi -infinite 
cylinder occupying X2<0. The first step in the solution of the problem is, 
therefore, the determination of the elastic properties E 2 , v 2 an ("infin- 
itely" long) equivalent shell of thickness h in terms of E-| , v-| , E' , v' , 
h, h 2 and R^ where E' and v‘ may or may not be the same as E-j and v-|. This 
step is briefly described in Appendix B where the composite cylinder is 
assumed to consist of fully-bonded thick-walled cylinders under axi symmetric 
plane strain condition and the equivalency of the radial displacements and 
axial strains is used to determine E 2 and v 2 . 

Details of the formulation of a homogeneous shallow shell containing a 
through crack by using a Reissner type transverse shear theory may be found, 
for example, in [9]- [11] and will not be repeated here. Thus, referring to 
Appendix A for the definition of normalized and dimensionless quantities, 
in terms of the displacement w, stress function <j> and the auxiliary functions 
ip and Q, the basic equations for the nonhomogeneous cylindrical shell shown 


in Fig. lb may be expressed as follows: 

vV(x i /u 1 .) 2 0= 0 , (i=l, y<0; i=2, y>0) (1) 

V 4 w+X?yf(l- Ki v2) 1 ^ = 0 , (i=l, y<0; i=2, y>0) , (2) 

= 0, (i=l, y<0; i=2, y>0) , (3) 

k. (1— v. ) 

— - o- 1 - V 2 n-n = 0, ( i =1 , y<0; i=2, y>0) (4) 
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where the constants , y- and k^. (i=l,2) are defined in Appendix A, v. is 
the Poisson's ratio, it is assumed that the cylindrical surfaces of the com- 
posite shell are free of tractions and the crack surface stress and moment 
resultants are the only nonzero external loads. The functions \p and n are 
related to the rotations as follows: 

_ d± + 12. „ _ _ K i^~ v i^ an 

K x " ax 2 ay * p y ay 2 ax * 

(i=l, y< 0; i=2, y>0) . (5) 


The normalized membrane, moment and transverse shear resultants are given by 

a 2 <j> >. _ a 2 4> 


N .... = , N„„ = 3 , N 


N xx ay 7 " 5 yy ax 
36 


xy 


ax ay 5 


M. 


xx 




= + . d-1. y<o-. 1 - 2 . y>o) , 


V = M + R V =~+R 

x ax B x * y ay p y 


aw 


(7) 

(8) 


Eliminating <j>, from (1) and (2) it may be shown that 
v 4 v 4 w + xjd-ic.7 2 ) |^T = 0 , ( i =1 , y<0; i=2, y>0) . 


(9) 


If we express the solution of (9) by 


w(x,y) = <*, 


J_ 

2tt 

2 

IT 


f 1 (x,a)e" lcty da + ^ f 2 (y,6)cos6xd6 , (y<0) 


2 f 


f-(y,6)cos6xd6 , (y>0) , 


(10) 


and assume the solution of the ordinary differential equations resulting 
from (9) and (10) of the form 
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f 2 (y»&) = R 2^) e 


OD 


m,x 

f-j{x,a) - R-j (a)e 

m r>y 

f 3 (y*B) = R 3 (3)e J 


m 2 y 


> 


The characteristic equations giving m^, m 2 and m 3 may be obtained as follows: 
m^-4a 2 m^+6a l+ m^-(4a 2 +K:-j A^Ja^m-j+a^Ca^+X^fl+K-ia 2 )] = 0 , (12) 

m 2 -(4s 2 +i<i A^)m|+(63 4 +K^ A^g 2 +x^)m 2 “4e 6 m 2 +B 8 = 0 , (13) 

m 3 “(4B 2 +K 2 A 2 )m 3 +(63 4 +tc 2 A 2 3 2 +A 2 )m2-43 6 m2 + ® 8 = 0 (14) 


We designate the roots of (1 2)-( 14) by m. • , i=l,2,3, j=l,...,8 and note that, 

• J 

properly ordered, they have the following properties: 


Re(m i j)<0, m i j +4 =-m. j , (i =1 ,2,3, j=l,2,3,4) . (15) 

Assuming that the composite shell shown in Fig. lb is loaded syrrmetrically, 
it is sufficient to consider the problem for x>0 only. Also, since the 
external loads acting on the crack surfaces are statically self-equilibrating, 
the functions f ^ , f 2 and f 3 must vanish at infinity and may, therefore be 
expressed as 


f-j (x,a) 

f 2 ( y >3) 

f 3 (y.B) 


E R,.(a)e 
j=l 13 


m lj x 


8 m~ .y 

E R ? ,(3)e ° 
j =5 ^ 

4 m-.y 

E R 3 ,(3)e 33 
j=l 33 


, x>0 , 

, y<o , 
, y>o , 


(16) 

(17) 

(18) 


Similarly, if we let 
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<J>(x,y) = i 


r X 

’ 

2u 

, 


OO 

! 1 

• 

; 7 r 



g-i (x,a)e lay da + — 

I TT 


g 3 (y, 3 )cosexde , (y>o) 


g 2 (y 5 e)cosBxdB, (y<0) 


(19) 


from (1), (2), (10) and (1 6)— (18) it may be shown that [10] 
X-. 4 R-i 4 ( a ) m-. . x 

g, (x,ot) = - 4- ct 2 z -4- — e lj , (x>0) , 

1 y l j=l p lj 


(20) 


A l 2 2 m 2j R 2i^) _ m 2.i y 


g 2 (y »3) = t i? = -fe 


1 j=5 p 2j 

2 


, (y<o) , 


xS 4 nio*Ro-i(s) 

g 3 (y»B) =4 I e 3j , (y>0) , 

J Mo j _ i 


2 j=l 


3j 


( 21 ) 


( 22 ) 


where 


p 1j = • p 2j = m 2r 62 • P 3j * 4 J' 82 


(23) 


Assuming now the solution of (3) and (4) of the form 


i ^ 7r 

n(x,y) = < 


h-j (x,a)e" 1ay da + f h 2 (y,3)sinBxdB , (y<0) 


2 
t r 


h 3 (y,B)sinSxdB , (y>0) 


(24) 


«l»(x,y) = < 


f J_ 

2tt 

2 

IT 


k-j (x,a)e~ lay da +| f k 2 (y,B)cosBxdB , (y<0) 


k 3 (y>B)cosBxdB , (y>0) 


( 25 ) 


It can be shown that 
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» 


h-j ( x ,a ) 

- A^aje , r-j -[a + » ( x> °)> 

(26) 

h 2 (y,e) 

= A 2 (e)e 2 , r 2 - [B 2 + (y<0) , 

(27) 

h 3 (y,e) 

= A 3 (s)e 3 , r 3 = -Ce 2 + (y>o) , 

(28) 

k-,(x,a) 

4 R-, £ (a) m, .x 
= s - e J » (x> 0 ), 

j=l K l p lj 1 

(29) 

k 2 (y,e) 

8 R ?i (0) nu-y 

= 1 J o -1 e d3 » (y<°) 7 

j=5 K l p 2j 1 

(30) 

k 3 (y,B) 

4 R-s (3) rn..y 
= S 7 - 5 — rr e 3J , (y>0) . 
j=l K 2 p 3j 1 

(31) 


The preceding formulation contains fifteen unknown functions R-^, 

( j-1 . »4) , ^2 j ’ 5 » • • • > 8 ) } R^ j , (j-1 . . ,4) , and , (i - l)2,3), which 

are determined from the following boundary and continuity conditions 


(see Fig. lb and Appendix A): 

u(x,-0)=u(x,+0), v(x,-0)=v(x,+0) , w(x,-0)=w(x,+0), (x>0) , (32) 

0 x (x,-O)=B x (x,+O), B y (x,-O)=0 y (x,+O), (x>0) , (33) 

E lN yy (x,-0)=E 2 N yy (x,+0), E 1 N xy (x,-0)=E 2 N xy (x,+0), (x>0) , (34) 

E 1 M yy (x,-0)=E 2 M yy (x,+0), E 1 M xy (x,-0)=E 2 M xy (x,+0), (x>0) , (35) 

B 1 V y (x,-0) = B 2 V y (x,+0) , (36) 

N xy (0,y)=0, M xy (0,y)=0, V x (0,y)=0, (y<0) , (37) 

N xx (0,y) = F-, (y) , (-d-j^-b-,) , 

u(0,y) = 0, (-co<y<_d 1 , -b 1 <y<0) , (38a, b) 
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(39a, b) 


M xx (0,y) = F 2 (y), ( — d-j <y<— b-j ) , 

3 x (0,y) = 0, (— <y<-dp -b 1 <y<0) , 

where and F 2 are the crack surface tractions obtained from the solu- 
tion of the uncracked shell under the given external loads. Note that 
for y>0 the assumed solution has the proper symmetry and gives 

N X y( 0 ,y)=Mxy(°,y)= v x (0 »y ) =3x (0 ,y)= u ( 0 ,y)= 0 , (y>o) . (40) 


It is seen that once the functions w, <}>, ip and ft are determined, e . , N. 

* * J 

. and , (i,j=x,y) may be expressed in terms of R - ^ and l\- by using 

(5)-(8). To complete the formulation of the problem the displacements u 
and v need to be determined. This may be done by using the Hooke's law 
and the following kinematic relations 


: ij " 2 (u i,j +U j,i +Z ,i U 3,j +Z ,j u 3,i ) ’ (l,J " 1,2) ’ 


(41) 


where the function Z(x-j,x 2 ) describes the middle surface of the shell. 

For the cylindrical shell under consideration Z 2 =0, Z -j -j =-l / R and refer- 
ring to Appendix A we find 



3N w 3N 

2(l+v. ) -M- 

v i ; 3y 3x 


3N 


+ v. 


XX 


9X 


X i 3 2 w 
vT 7 ay 7 


(42) 


3V 


= N...-V.N 


ay yy i xx * 

where i=l for y<0 and i=2 for y>0. From (42), (43) and the formulation 
given in this section it can be shown that 


(43) 


1 X 1 

u(x,y) - 


4 (2+v, )a 2 -m, . m,.x-iay. 

y 1 .?3 D f ra \ m p 'J 


j=l 




2 x l 2 r°° 8 (2+v.)m^.-3 2 m ? .y 

1 i - - 7 -=^ R 9 .(B)3e J sinexdp 

0 J=5 p 2j ^ 


+ - — 7 
^ j 


A 1 

+ — 2 : x w(x,y) , (x>0, y<0) , 
y l 


(44) 
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9 \ 9 2 r 4 (2+v ? )m 2 .-B 2 nu.y 

u(x,y) e 52-^ Ro n -(e)ee ^ sinBxdB 

11 y 2 J q j=l p 3j 33 

X 2 Z 


+ 77 x w(x,y) , (x>0, y<0) , 
y 2 


»«*•») --srS 


2 « 


m 1 .x-iay 


E -—-2 H. •(a)(m?.+v 1 a 2 )e 
3=1 P lj 3 13 

n Xi ^ r 0 . nu, (Ha *y 

■777 E 33 cossxdB, 

77 y l J j=5 p 2j 33 1 33 


v(x,y) = — 

it y 2 


2 *2 f 4 rn 


(x>0, y<0) , 

m ^-y 


4 "lo,- '"o-jj 

E # R^(8)(B 2 +v 9 m 2 .)e J cosBxdB , 
0 j=l P 3j 33 1 33 

(x>0, y>0) . 


( 45 ) 


(46) 


(47) 


3. The Integral Equations 

Referring to the general formulation of the problem given in the 
previous section it is seen that the first thirteen conditions (33) -(37) 
are homogeneous and may be used to eliminate thirteen of the fifteen 
unknown functions and ( i =1 ,2,3; j=l,...,4; k=l,2,3).The mixed boun- 
dary conditions (38) and (39) would then give the integral equations to 
determine the remaining two. To derive the integral equations we first 
introduce the following new unknown functions which are the complements of 
the known crack surface tractions F-j and F 2 : 

G-j (y) = -^-u(+0,y) , G 2 (y) = e x (+0,y) , (y<0) . (48) 


From the solution given in the previous section and from (48) it may be 
shown that 


Sty) - - i 


tt y. 


5 (H 4 ) ° 2 -P 1 .i 


I 

3=1 


Plj 


Rlj(a)ae la ^da , 


(49) 
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G 2 (y) = 


i 

2 TT 


4 am 1 .R-j • (a) 
z — - n -- -i e lay da 

j-l K l p lj-’ 


tc-j(l-v-j) f 


4 TT 


a 2 A-| (a)e lay da 


( 50 ) 


By inverting (49) and (50) and by using (33 )- (37 ) one may express all fif- 
teen unknown functions and in terms of G-j and G 2 . From (38b), 

(39b) and (48) it may be observed that 


G 1 (y) = 0, (i=l ,2; -<*><y<-d-, , -b 1 <y<0) 


(51) 


and (38b) and (39b) would be identically satisfied if 


r b i 


j i 


G^yjdy = 0, ( i =1 , 2 ) . 


(52) 


From the formulation of the problem the conditions (38a) and (39a) may be 
expressed as 


N (+0,y) = lim (~) I ^z~ >'i j 

xx x-vfO ^ U 1 J j=l p lj IJ 


4 4 m,.x-iay 

r ^IW=Oe lj 


da 


p 2 

+ l (-1 

7 T V U -| 


4 

8 m oA '"2i J 

z Rp,-(B)e J cosBxdB = F, (y) , 
j=5 P 2j ^ 


m,,y 


( _ d-|<y<-b-j ) 


(53) 


M xx^ +0,y ^ = llm ^27 
xx x->+0 y l c 


4 mfj-v^ 2 


tt^t R ij (a)e 
j=l l p lj 1 IJ 


m l j x " 1 ay 


iici(l-vi )' 

4tt 

o 

<1 O-vi )‘ 


r,x-iay or a p -vomp. 1 

a^A^cJe 1 fa-- j [ r. ■ r ] -^j(B) e 


da 


f r 5 B2 : ^2i 

O J=s 1 ' 2j 


m 2j y 


r 2 y 


Br 2 A 2 (e)e jcosBxdB) = F g (y ) , (-d-j <y<-b-j ) . (54) 
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If we now write R. . and A. in terms of 6-j and G2 and take into 

J 

account (51), (53) and (54) become 
-b. 


'1 


-d. 


i H k j(y,t)Gj(t)dt = F k (y), (k=l,2, -d,<y<-b.|) , 
d 1 


( 55 ) 


H ki (y,t) = lim { 
KJ x-»+0 


B^U.oOe 1 ^ t-y ^ a dct + 


C k j(y,B)cosBxdB) , (56) 


where B kj . and C k j, (k,j=l,2) are known functions. The derivation of these 
functions are rather lengthy but quite straightforward and will not be 
reproduced in this paper. From the viewpoint of obtaining the correct 
singular behavior of G-j and G£ near and at the end points y=-d^ and y=-b-j 
and a sufficiently accurate solution of the system of integral equations 
(55), it is essential that the dominant parts of the kernels given by (56) 
be separated. This can be done by examining the asymptotic behavior of 
B k j(x,a) and C kJ -(y,B), (k,j=l,2), for }a|-*» and B-*», respectively, by 
separating the asymptotic terms Bj. and C k j m and by evaluating the corre- 
sponding integrals in closed form. A key step in this process is the deter- 
mination of the roots of the characteristic equations (12)-(14) (see, 
for example, [10]). It may then be shown that for large values of |a| and 
8 the roots have the form 

2 


m 


Id 


m 2j = ed 


. la 1(1 + - dj- + ...). 


2a' 


+ !2j M + 
+ w W " 


) , (d=5 8) 


m 3j “ e(1 + 2B 


3d 


P3d 

W 


•••) > (d = l»*»*>4) 


(57) 

(58) 

(59) 


Also, from (26)- (28) it may be seen that for |a|»l and B>>1 we have 
^i = -I a l(l + — ), r 2 = 8(1 + (1-v 1 ) " **•) ’ 

r 3 = “ B(1 + ^TT^TF " * 


(60a-c) 
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After evaluating the dominant kernels the integral equations (55) 
may be expressed as 
-b. 


'1 


{[k ls (y,t)+k 11 (y J t)]G 1 (t)+k 12 (y,t)G 2 (t)}dt = 2irF 1 (y), (-d 1 <y<-b 1 ) 


-b. 


(1-v-, 2 ) {k 21 (y,t)G 1 (t)+[k 2s (y,t)+k 22 (y,t)]G 2 (t)}dt 


~dn 


hy 


"Ml 

= 2lr a F 2 (y)» (‘d^y^b-j) , 


(61) 


(62) 


where 


k ls (y, t )=k 2s (y,t)= Jp- 


- ( 


c,(3c 3 +2 c 4 )+2 c 2 c 5 , 6c 


c l c 6 


.) _ J — + (_jL) 
' t+y v c r ‘ 


Tt+yF 


( 4c 3, t y 2 
c c (t+y) 3 * 


(63) 


c-|=-3+v 2 -(1+Vi )E 2 /E-| , c 2 =1+E 2 /Ei , C2 = l+v 2 -(l+v-| )E 2 /E-j 


c 4 =-v 2 +ViE 2 /Ei , Cg-3-v 2 - (3-v-j ) Eg/ E-j > Cg 1 +v 2 +(3-v-j ) E 2 /E-j 


(64) 


and the kernels k-.(y,t), (i,j=l,2) are bounded in the closed interval 

* vJ 

- d il(y,t)<-b 1 (including the case b-j =0) . The dominant or singular kernels 
k.j s and the bounded kernels k^., (i,j=l,2) are obtained from the integrals 
of the form: 


k ls (y,t) 


lim { 
X-++0 


B llm (x s a)e l(t - y)a da + 


C-| lco (y,B)cos6xd6}, (65) 


— CO 
oo 

+ j [C lj -(y,6)-C lj . oo (y,B)]dB , ( j=l ,2) . (66) 

o 
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The important point to observe about the integral equations is that, as in 
the homogeneous shells, only the diagonal kernels contain singular terms 
and the dominant kernels k^ s and k^ s corresponding to membrane and bending 
loads are identical. This physically expected result is possible again 
because of the use of a transverse shear theory in formulating the problem. 
Also, it can easily be shown that the dominant kernel k^ s found for the 
nonhomogen eous shell in this study is identical to that obtained for two 
dissimilar bonded half planes with a crack perpendicular to the interface 
under plane stress conditions which is given in [14]. The plane stress 
problem is, of course, a limiting case of the shell problem and this, too, 
is the expected result. 


4. The Stress Intensity Factor 


After solving the integral equations, clearly any desired field quan- 
tity may be obtained from integrals with appropriate kernels having G-j 
and G 2 as the density functions. For example, it may be observed that 
before going to the limit, (53) and (54) give N (x,y) and M (x,y) every- 
where in the shell. In particular we note that (55) or (61) and (62) are 
valid for x=0 outside as well as within the cut -d-|<y<-b.j. Thus, through 
a simple asymptotic analysis of (61) and (62) one can obtain the stress 
intensity factors at the crack tips which, for b>0, are defined by (Fig. 
lb) 


k. (-b,Xo) = lim 
1 J x 2 ->-b+0 


i/2(x 2 +b) a 11 (0,x 2 ,x 3 ) , 


(67) 


k-i(-d,Xo) = lim /-2(x<j+d) a, ■, (0,x o ,x,) . (68) 

1 J x 2 ->-d-0 6 11 * J 

In the shells the in-plane stress components are obtained by combining 
membrane and bending stresses as follows: 


b 

°i j( x l » x 2* x 3^ = CT ij + °i j 

m 1 M % b 

a ij " TT N ij( x i > x 2^ » a ij 


, (i »J = 1 »2) 


12x 3 

~W~ M ij^ x l ,x 2^ 


(i,j=l,2) 


(69) 

(70) 
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Now for b>0 from (61)- (63) it is seen that the dominant kernels k,_ and 

“1 ^ 
k 2s consist of (t-y) only and hence the solution of the integral equations 

may be expressed as [15] 


Gj(y) 


p .lM 

[-(y+dpfy+b,)^ 


(j=l,2; -d-j <y<-b-| ) 


(71) 


where P-j and P 2 are unknown bounded functions. By substituting from 
Appendix A and (71) into (61) and (62) and using (67)-(70), it can be 
shown that 


k-j (-b,x 3 ) = 


‘I 


- ^-lim /-2(x ? +b) — [u, (+0,x ? )+x^B 11 (+0,x ? )] 
d x 2 +-b-0 d 9x 2 1 c 6 11 c 


Y /a [P, (-b, ) +y- p 2 ( ' b l );1 ’ 


(72) 


k-,(-d,x 3 ) = 


-J- lim 

x 2 ->-d+0 


72(x 2 +d) [^1 (+0 ,x 2 )+x 3 3-j -j (+0 ,x 2 ) ] 


■ -f & [P n (-d, ) + ^ p 2 (- d i)] 


(73) 


In the case of b=0 from (63) it is seen that the dominant part of 
the kernel is a generalized Cauchy kernel, that is it contains, in addi- 
tion to (t-y)" 1 , terms which become unbounded as the variables y and t 
approach the end point -b-|=0 (Fig. lb). The contribution of these terms 
to the singular behavior of the solution at y=0 can be studied by assuming 
the solution of the integral equations (61) and (62) as 


Gj(y) - Sj(y)(-y) Y (y+d-j )“ , (-l<Re(y,w)<0) , (74) 

and by following the function theoretic method (see, for example, [14] 
and [15]). Thus, by substituting from (74) into (61) and (62), the charac- 
teristic equations giving y and w may be obtained as follows: 

2c~ 3c,+2c>, 2 c 9 Cc 

cos-jry - — — y(y+2) - ( — + --■■■ - - -) = 0 , (75) 

c 6 c 6 c l c 6 
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COSiru = 0 


(76) 


At the crack tip y=-d-| which is embedded in a homogeneous medium (76) 
gives u=-l/2. The characteristic equation (75) found for y is identical 
to that given in [14] for the plane stress case and its examination would 
show that for all material combinations the equation has only one root 
satisfying -l<Re(y)<0 and this root is always real. In [14] it was found 
that for the plane problem in the small neighborhood of the singular point 
y=0 located at the interface the stress state has the form 


k. 

a,-,- (r,e) = — r Y f - - ( 0 ) , i,j=r,e) , 
ij P) ij 


(77) 

I 


where r and e are the polar coordinates and the functions f . . are dependent 

* 

on the bimaterial constants E^, v k , (k=l,2) and are given in [14]. The 
constant k-| is again defined as the stress intensity factor and is obtained 
from the calculated values of a_ Q (r,0) by normalizing f Qn (0)=l. Thus, 

DU U U 

in the present shell problem the stress intensity factor at the crack tip 
y=0 may be defined as 


k, (0,Xo) = lim 
1 J x 2 -^0 


'fl X2’" Y a-|-|(0,x 2 ,x 3 ) . 


(78) 


To evaluate the stress intensity factor k-](0,x 3 ) the asymptotic expression 
of cr-j i for x 2 >0 is needed. This can again be obtained in terms of G-j and 
G 2 from the basic formulation of the shell given in this paper. After 
somewhat lengthy but straightforward analysis it may be shown that 


Eo 3c f -+c 1 +2y(ci+C/-) „ 

k (0,x 3 ) = - { - 6 -1 ■ J: 6 ) lim (-x 2 )' y 

6 JZ ^CgSinuy x ^+0 1 3x 2 


[u 1 (+0,x 2 ) 


x 3^1 1 ( + ® ,x 2^ 


E 2 a" Y ( 


3 c 6 + c i+ 2 ( c i+c 6 )y 

C-jCgSinny 


KS^(O) +-fs z m 


(79) 
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Also, at the crack tip x 2 =-d the expression 
E 1 _ a 

k 1 (-d,x 3 ) = -j- lim > / 2(x 2 +d} — [u-,(+0,x 2 ) + x 3 b 11 (+0,x 2 )] (80) 

i ■' x 2 -+-d+0 *2 " ' 

is still valid and by using (48) and (71) can be written as 


k-|(-d,x 3 ) & [Si (-d-j ) + a 3 s 2 (-d-j ) H . 


(81) 


5. Results 

For b>0 the integral equations (61) and (62) subject to (52) may 
easily be solved by using the Gauss-Chebyshev quadrature formulas described 
in, for example, [14] by assuming the solution in the form (71). The 
stress intensity factors may then be obtained from (72) and (73). From 
(72) and (73) it is seen that the values of the bounded functions P-j and 
?2 a t the end points -d-j and -b-| are associated with respectively the 
membrane and the bending components of the stresses near the crack tips. 

For the pressurized shell shown in Fig. lb, one can, therefore, define the 
following normalized stress intensity factors: 


k_(r .) 
m j' 


Mr-,0) k. (r- ,h/2)-k, (r- ,0) 

— ! — J } k ( r .) = — J ! — J 

(pP^/hJfa" b J (pR^/h) Ja 


(82) 


where r-j=-b and r 2 =-d and the stress intensity factors k^(r^,x 3 ), (j=l ,2) 
are given by (72) and (73). For b=0 the definition of the normalized 
stress intensity factors at r 2 =-d would remain the same as in (82) and 
may be obtained from (81). At r^=-b=0 we have 


M0,0) k, (0,h/2)-k 1 (0,0) 

MO = — 1 , k b (0) = 1 1 

m \ "Y d 


(pR.j/h)a 


-Y 


(pRj/h)a 


-Y 


(83) 


where k-j is given by (79). In (82) and (83) R-j is the inner radius of the cylinder. 

The reinforced shell shown in Fig. la is considered as an example. 
The main shell having properties ,v^ is steel and the reinforcing shell 
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is a fiber reinforced composite. The dimensions of four different 
cylinders used in the analysis are shown in Table 1. The table also 
shows the elastic constants E 2 ,v 2 the equivalent homogeneous shell 
shown in Fig. lb. The derivation of E 2 ,v 2 anc * the so ^ ut ^ on the nonhomo- 
geneous shell in the absence of any cracks giving the crack surface 
loads F-| and F 2 are given in Appendix B (see, (61) and (62)). It should 
be noted that F^ and F 2 are functions of y=x 2 /a. One may also note that 
the effective elastic constants E 2 ,v 2 as well as F-j,F 2 are dependent on 
the axis! constraint in the cylinder. In the examples it is assumed that 
either there is no axial constraint (labeled as a cylinder with "open 
ends") or the ends of the cylinder are "closed", corresponding to the 

2 

total axial force in the internally pressurized cylinder P=0 and P=nR-p, 
respectively (Eq. B6). 


Table 1. Dimensions and the effective material constants 
(E 2 ,v 2 ) of the composite shell used in numerical 

examples. 


Shell # 


1 

2 

3 

4 

E ( psi 

t l *GPa 

V 1 

E 1 (reinforcing shell){^ 

v' (reinforcing shell) 
h (thickness) { f | i n * 

R (mean radius) -C^ 0, 

h 2 (reinforcing shell)!™* 
R/h 

3x1 0 / 
207 

0.3 

2xl0 7 

138 

I 0.1 
! 1 
0.0254 

100 

2.54 

1 

0.0254 

100 

3xl0 7 

207 

0.3 

2x10 7 
138 

0.1 

0.615 

0.0156 

23.6925 

0.6018 

0.615 

0.0156 

38.524 

3x1 0 7 
207 

0.3 

2xl0 7 

138 

0.1 

0.5 

0.0127 

17.75 

0.4509 

0.5 
1 0.0127 

35.5 

~ 3x1 0 7 
207 

0.3 

2xl0 7 

138 

0.1 

1 0.404 

1 0.0103 

11.798 

0.2997 

i 0.404 

l 0.0103 

• 

29.203 

Open Ends: 

c 

rPSi 

l GPa 

50.2xl0 6 

49.72xl0 6 

49.66xl0 6 

49.48x10 s 

b 2 

346 

343 

342 

341 

v 2 


0.2214 

0.2174 

0.2168 

0.2153 

-Y 


0.449518 

0.450677 

0.450827 

0.451265 

Closed Ends: 

r 

rPSi 

l GPa 

50.06x10 s 

49.37xl0 6 

49.28x10 s 

49.03x10 s 

E 2 

345 

340 

340 

338 

v 2 


0.2272 

0.2321 

0.2327 

0.2345 

Zl 


0.440538 

0.450728 

0.450886 

0.451323 
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For these two cases the values of E 2 and v 2 calculated from Appendix B are 
given in Table 1. 

The stress intensity factors obtained by using dimensions and proper- 
ties of four different shells shown in Table 1 are given in Tables 2-5. 

The normalized stress intensity factors given in the tables are defined 
by (82) for b>0 (or c>a) and by (81) and (83) for b=0 (or c=a). One may 
observe that as the relative distance of the crack to the interface c/a 
increases the stress intensity factors approach those given for the homo- 
geneous shell [9]. The tables also show that the membrane component k m 
of the stress intensity factor (which is by far the dominant part) decreases 
with the decreasing crack distance to the boundary. Even though the 
results found for "open" and "closed" ended cylinders are different, for 
the axial through crack geometry under consideration the differences seem 
to be relatively insignificant. It should be emphasized that the tables 
show the stress intensity factors normalized with respect to (pR^/hJ/a" 
which is the corresponding flat plate (or plane stress) result under the 
same membrane loading as the shell. Thus the variation in the stress 
intensity factors as a function of a/h and R/h is entirely due to curvature- 
and thickness effects. 

In Tables 2-5 the membrane and bending components of the stress inten- 
sity factors are given separately. In all cases, the calculated bending 
components were such that the stress intensity factors on the outside sur- 
face x 3 =h/2 were greater than that on the inside surface. Figures 2-5 
show the stress intensity factor k^(c 1 -,h/2), (c..=-d,-b,Q) obtained from 
(72), (73), (79) and (81) for some selected shell-crack geometries. The 
effect of reinforcement may be clearly observed from Figures 2 and 3. The 
asymptotic behavior of k-j (-b ,h/2 ) as c+a is due to different definitions of 
k-j for b>0 and b=0 as given by (67) and (78). For these two cases the 
cleavage stress in the close neighborhood of the crack tip may be expressed 
as (Fig. lb) 

M-b,h/2) 

a 11 (0,x 2 ) = , (r=x 2 +b, b>0, x 2 >-b) , (84) 

/2r 
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> 


(85) 


. k-i (0,h/2) 

^n(0»Xo) = ~ 

11 2 /2 r Y 


k 1 (0,h/2)r %+Y 

y/Zr 


(r=x 2 >0, b=0). 


Thus, since l/2+y>0 (see Table 1), the stress intensity factor defined on 
the basis of the conventional square root singularity becomes 


lim k.(-b, h/2) = lim a-,-, (0,x 9 ) = lim k,(0,h/2)r z Y = 0 . (86) 

b^O r->0 11 6 r->0 1 

b+O 

Figures 4 and 5 show the normalized stress intensity factors on the 
outside surface Xg=h/2 for b=0 obtained from (79) and (81). Here, too, 
it may be seen that even for a relatively very shallow shell (shell no. 1), 
there is considerable curvature and thickness effect. 
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Table 2. Membrane component k m of the normalized stress 

intensity factor in a pressurized composite shell 
with open ends which contains an axial crack. 


Shell 


k (-b) 
m ' 

k(-d) 
m ' 

No. 

c /a 

a/h 

1.0 

1.1 

1.5 

2 

10 

1.0 

1.1 

1.5 

2 

10 


1 

1.812 

H 

Qgj?| 

0.809 

m 

0.789 

0.798 

0.820 

0.839 


1 


1.878 

o 


0.864 

1.050 


0.852 

0.885 

0.916 

1.049 


1.973 

0.759 

0.874 

0.934 

1.082 

0.903 

0.918 

0.960 

0.998 

1.079 


10 

3.097 

1.203 

li 447 

1.534 

1.542 

1.442 

1.472 

1.532 

1.549 

1.542 


1 

1.849 

0.721 


EgSI 

1.041 

0.830 

0.840 

0.869 

0.895 

1.041 

2 


1.998 

0.782 

1 

BB 


0.929 

0.944 

0.988 

1.027 

1.097 



2.198 

0.862 

1.010 

in 

1.175 

1.042 

1.061 

1.113 

1.152 

1.176 


10 

4.426 

1.730 

2.057 

2.101 

2.096 

1.999 

2.036 

2.091 

2.096 

2.096 


1 

1.853 

0.724 

0.817 

EH 


0.835 

0.845 

0.875 

0.902 

1.049 

3 


2.013 

0.788 

0.910 

Ha 

1.106 

0.940 

0.956 


1.039 

1.104 



2.227 

0.876 

1.027 

1.107 

1.189 

1.058 

1.078 

1.131 

1.170 

1.189 


10 

4.581 

1.790 

2.125 

2.164 

2.161 

2.062 

2.100 

2.155 

2.160 

2.161 


H; 

1.866 

Esa 

0.829 

0.871 

1.057 

0.848 

0.860 

0.890 

0.919 

1.058 



2.053 

M 

0.937 

1.004 

1.124 

0.969 

0.985 

1.031 

1.07V 

1.123 

4 

1 

2.304 


1.073 

1.156 

1.225 

1.102 

1.123 

1.177 

1.214 

1.226 


10 

4.980 

1.947 

2.297 

2.327 

2.326 

2.225 

2.264 

2.319 

2.324 

2.326 
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Table 3. Bending component k^ of the normalized stress intensity 
factor in a pressurized composite shell with open ends 
which contains an axial crack. 


Shell 


k m<- b > 

k(-d) 

m 

No. 

c/a 

a/h 

1.0 

1.1 

1.5 

2 

10 

1.0 

1.1 

1.5 

2 

10 


1 

0.005 

0.002 

0.008 

0.014 

0.039 

0.015 

0.016 

0.020 

0.024 

0.037 

1 

2 

0.047 

0.021 

0.035 

0.044 

0.046 

0.042 

0.044 

0.050 

0.055 

0.043 

3 

0.087 

0.042 

0.062 

0.075 

0.064 


0.072 

0.078 

0.082 

0.064 


10 

0.225 

0.134 

0.202 

0.212 

0.197 

0.191 

0.196 

0.204 

0.203 

0.198 


mm 

0.033 

0.015 

0.026 

0.034 



0.035 

0.041 

EB 

0.039 

2 


0.101 


0.069 

0.081 


Ell" 

0.079 

0.085 

EB 

0.070 


1 

0.154 

0.074 

0.108 

0.123 

0.110 

0.111 

0.114 

0.121 

EB 

0.110 


10 

0.164 

0.124 

0.211 



0.207 

0.198 

0.200 

0.202 

0.207 

0.203 

0.198 


1 

0.037 

0.017 



0.044 

0.036 

0.038 

0.043 

0.047 

0.040 

. 3 


0.107 

0.050 

ESI! 


0.074 

0.079 

0.082 

0.089 

0.092 

0.074 



0.161 

0.078 

0.117 

0.128 

0.114 

0.115 

0.119 

0.125 

0.127 

0.114 


10 

0.148 

0.108 

0.197 

0.200 


0.194 

0.197 

0.201 

0.197 



0.193 


1 

0.046 

0.021 

0.032 

Ksa 


0.041 

0.043 

0.049 


0.042 



0.120 

0.055 

0.082 

M 


0.088 

0.091 

0.097 

0.099 

0.082 

4 


0.176 

0.086 

0.124 

Em 


0.126 

0.129 

0.135 

0.136 

0.125 


10 

0.102 

0.099 

0.186 

0.172 

0.172 

0.178 



0.177 

0.172 
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Table 4. Membrane component k m of the normalzied stress intensity factor 
in a pressurized composite shell with closed ends which contains 
an axial crack. 


Shell 


k (-b) 
nr ' 

! k m<- d > 

No. 

c/a 

a/c 

1.0 

1.1 

1.5 

2 

10 

1.0 

1.1 

1.5 

2 

10 


1 

1.909 

1 

0.820 

0.848 

0.985 

0.828 

0.836 

SI 

0.874 

0.994 

1 


1.969 


0.856 

0.897 

1.047 

0.872 

0.883 

SB 

0.938 

1.047 



2.059 

0.795 

0.907 

0.959 

1.081 

0.929 

0.942 

SB 

1.011 

1.078 


10 

3.158 

1.229 

1.458 

1.533 

1.542 

1.448 

1.475 

1.531 

1.546 

1.542 


1 

1.932 

0.760 


0.887 

1.039 

0.865 

0.874 

0.900 

0.923 

1.042 

2 


2.074 

0.818 

0.932 

0.986 

1.099 

0.954 

0.968 

1.006 

1.039 

1.097 



2.265 

0.897 

1.037 

1.104 

1.176 

1.059 

1.077 

1.123 

1.156 

1.176 


10 

4.446 

1.750 

2.059 

2.098 

2.096 

2.001 

2.036 

2.089 

2.095 

2.096 


H 

1.936 

0.764 

0.856 

BH 

1.044 


0.879 

0.905 

0.929 

1.046 



2.087 

0.825 

0.942 

US 

1.106 

0.964 

0.978 

1.018 

1.051 

1.104 

3 

1 ‘S' 

2.293 

0.910 

1.053 

1.122 

1.189 

1.075 

1.093 

1.139 

1.173 

1.189 


10 

4.593 

1.810 

2.126 

2.161 

2.161 

2.063 

2.100 

2.153 

2.159 

2.161 


1 

1.944 

0.773 

0.867 

0.905 

1.055 

0.882 

0.892 

0.920 

0.944 

1.056 

4 


2.122 

0.846 

0.968 

1.026 

1.124 

0.990 

1.006 

1.047 

1.080 

1.123 



2.363 

0.946 

1.098 

1.169 

1.226 

1.117 

1.136 

1.184 

1.216 

1.226 


10 

4.977 

1.966 

2.297 

2.324 

2.326 

2.227 

2.265 

2.318 

2.324 

2.326 
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Table 5. Bending component of the normalized stress intensity 
factor in a pressurized composite shell with closed ends 
which contains an axial crack. 


Shell 



• 

■k b (-b) 





k b (-d) 



No. 

c/a 

a/h 

1.0 

1.1 

1.5 

2 

10 

1.0 

i.i 

1.5 

2 

10 


1 

0.010 

0.003 

0.010 

0.014 

0.034 

0.015 

0.016 

0.019 

0.022 

0.032 

1 

2 

0.051 

0.023 

0.035 

mm 

0.045 

0.041 

0.042 

0.047 

0.051 

0.042 

3 

0.091 

0.042 

0.061 


0.061 

0.067 

0.069 


0.078 

0.064 


10 

0.230 

0.135 

0.199 

0.209 

0.197 

0.190 

0.195 

0.203 

0.203 

0.198 


1 

1 

0.016 

0.026 

0.033 

0.040 

0.033 

0.034 

0.038 

0.041 

0.037 




0.048 

0.068 

0.079 

0.071 

0.074 

0.076 

0.082 

0.085 

0.070 

2 


1 

0.076 

0.107 

0.120 

0.110 

0.109 

0.112 

0.118 

0.119 

0.110 


10 



0.171 

0.125 

0.219 

0.205 

0.198 

0.199 

0.202 

0.207 

0.203 

0.198 


mm 

0.042 

0.019 

0.027 

0.035 

0.041 

0.034 

0.036 

0.040 

0.044 

0.038 


Kfl 

0.111 

0.051 

0.072 

0.083 

0.074 

0.078 

0.079 

0.085 

0.088 

0.074 

3 

1 • 


0.080 

o.m 

0.125 

0.114 

0.114 

0.117 

0.122 

0.123 

0.114 


10 

B 



0.200 

0.192 

0.195 

0.197 


0.197 

0.193 


1 

0.050 

0.022 

0.033 

in 

0.045 

0.039 

0.042 

0.046 


0.041 



0.125 

0.058 

0.081 


0.083 

0.086 

0.089 

0.092 

0.096 

0.083 

4 



0.088 

0.123 

IB9 

0.126 

0.124 

0.127 

0.134 

0.134 

0.125 


10 

0.112 


0.184 

0.179 

0.172 

0.179 

0.180 

0.181 

0.177 

0.172 
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APPENDIX A 

Normalized quantities used in the formulation 
of the composite shell problem 


x = x-|/a, y = x 2 /a, z = x 3 /a, b-j = b/a, c-j = c/a, d-j = d/a , (Al) 

u = u-j/a, v = u 2 /a, w = u 3 /a, e x = = B 2 , (A2) 

r F(x-j ,x 2 )/(E 1 ha 2 ) , y < 0 , 

<f>(x,y) = < (A3) 

[ F(x 1 ,x 2 )/(E 2 ha 2 ) , y > 0 , 

f°i j/Ei » y < 0 » 

a (x,y,z) = < (a,6=x,y; i,j=l,2) , (A4) 

ho /E 2 ■ » * 0 • 

r N^/C^h) , y < 0 , 

N o(x,y ) =<, (a,8=X,y; i ,j=l ,2) , (A5) 

aB [ N ij /(E 2 h) , y > 0 , 

c M^/C^h 2 ) , y < 0, 

M R (x,y) =< (a,8=x,y; i,j=l,2) , (A6) 

l M ij/(E 2 h 2 ) , y > 0, 

r V k / (B-j h ) , y < 0 , 

V (x,y) = < (y=x,y; k=l,2) , (A7) 

Y ( V k /(B 2 h) , y > 0 , 

B i = 5E./(12(l+v i )) , ^ = E./(B iW J) , (i=l ,2) , (A8) 

u- = 12(l-v?)a 2 /h 2 , xj = 12(l-v?)a 4 /(R 2 h 2 ) , (i=l,2) . (A9) 
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APPENDIX B 

Solution of the uncracked composite shell problem 

The first problem here is the determination of elastic constants E 2 
and of a cylindrical shell which has the same stiffness as a long layered 
cylinder (see Fig. la and lb for x 2 »0). The second problem is the solu- 
tion of the bonded shells shown in Fig. lb in the absence of any cracks 
in order to determine the crack surface tractions F-j (y) = (0,x 2 )/(hE-j ) 

and F 2 (y) = M-j -j (0,'x 2 )/(h 2 E-| ) which are used in the perturbation solution 
of the cracked shell as the input functions. In the first problem we 
use the following basic (axi symmetric) solution of a thick-walled cylinder 
under plane strain conditions: 

a (r) = A + -~r , o ( r ) = A - , 

rr r z 90 r z 

w(r) = 1 ~ v ^ 2v Ar - ^ ^ - ve zz r , (a<r<b) , (Bla-c) 

where A and B are unknown constants and w is the radial displacement. 

For the three cylinders we have (Fig. 1) 

a-j = R.j , b^ = R^+h, A = A-j , B = B-| , \> = , E = E-j , (B2) 

a 2 = R i +h ’ b 2 = R 1 * +h+h 2’ ^ = A 2’ B = B 2* v = v *> B = E 1 > (B3) 

a^ = R. , b^ = R^+h, A = A^, B = B^, v = v 2> E = E 2 , (B4) 

where the dimensions , h, h 2 and the elastic constants E-j , v-j , E’ 

and v 1 are known. The bounday and continuity conditions 

a lrr^ a l ^ = “ p * a 2rr^ b 2^ = °» a lrr^ b l^ = °2rr^ a 2^ » 

w-, ( b-j ) = w 2 (a 2 ) , CT 3rr (a 3 ) = -P , o 3rr (b 3 ) = 0 (B5) 

are used to eliminate A. and B. , (i=l,2,3). Additional information 
needed to account for e zz is 
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'zz 


(a zz-' ,a rr'”e6 )/E ’ 


°zz dA 


= p 


(B6a,b) 


where P is the total axial force. If the cylinder is fully constrained 
e zz =0, if there is no axial constraint, P=0 and, if the ends of the 

cylinder are closed and if it is subjected to an internal pressure p, 

2 

then P = irR^p. In this problem it is assumed that the cylinders 1 and 
2 are fully bonded along r=b 1 =a 2 giving e i zz =e 2 zz * 

The two shells are said to be equivalent if 


£ 3zz e lzz 


w 3 (R) = w-j(R) , 


(B7a,b) 


where R = R^+h/2. After eliminating A. and (i=l,2,3), (B7a) and (B7b) 
would give the equivalent (or effective) elastic constants E 2 and v 2 in 
terms of , E', v‘, R^ , h and h 2 . 

The next problem is the evaluation of the stresses in the bonded 
shell shown in Fig. lb in the absence of any cracks. Referring to [16] 
and Fig. lb, the solution of the "open-ended" or axially unconstrained 
cylindrical shell problem may be expressed as 


u 3 (x 2 } = \ 


&1 *2 dR 2 
e (C-jCos6^x 2 +C 2 sinB-|X 2 ) - , x 2 <0 


®2 X 2 dK^ 

e (C^cosB^+C^sing^ - x 2 >0 , 


>R 2 


where 


6 1 = 3(l-v 2 )/R 2 h 2 , e 2 = 3(l-v|)/R 2 h 2 


(38) 


(B9) 


and the constants C-j,...,C^ are determined from 

u 3 (-°) = u 3 (+0) , u 3 (-0) = u 3 (+0) , (BIO) 

D 1 4 U 3 ( “° ) = D 2 ^ U 3 (+0) ’ °1 4 U 3 ( ' 0) = D 24 U 3 (+0) ’ (B11) 
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(B12) 


E^h 3 E 2 h 3 

D 1 = 12{l-vf) ’ °2 = 12(1 * 

Note that (BIT) represents the continuity of M 22 and V 2 . Equations 
(B8), (BIO) and (BIT) may be shown to reduce 


C r C 3 = P R 2 (E 2 -E 1 )/(hE 1 E 2 ) , 

C3 ) Qyl B-] , 

C 2 = -D 2 s 2 C 4 /( D 1 B 1 ) » 

C 2“ C l = (C 4 +C 3 )D 2 b|/(D 1 3 3 ). 

After determining C-|,..„,C^ from (B13), by 
E l h 

1 ( x 2 ) - “ ^ u 3 ( x 2 ) , M-j 1 ( x 2 ) - - v 


observing that [16] 

1 D 1 Ix^ u 3 (x 2 ) ’ 

(x 2 <0) , 


(B13a-d) 


(B14) 


the input functions F, and F 9 may be expressed as (see (38a), (39a) and 
Appendix A) v ' 


My) - i 


i vy 


[C-, cose, ay 


CjSing-j 


ay] 


hEi 


(y<o) 


F 2 (y) = 


2v 1 D 1 

w 


Si ay 

e [-C-j sing-jay + C 2 cosg-|ay] , (y<0) . 


(B15) 

(B16) 


Similarly, if the ends of the pressurized composite cylinder shown 
in Fig. lb are closed, following [17] we find 


^ ) Note that F-j and F 2 used in the perturbation problem, correspond to 
the membrane and bending resultants obtained from the uncracked 
cylinder with the opposite sign. 
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f Pi x o v l - ^ nR^ 

e [B-jCOS0.jX 2 +B 2 sine-jX 2 ] + (— g— ) , (x 2 <0) 

-BpXp v 2 “2 n D 2 

e [B^ co s 6 2 x 2 +B^s i n 0 2 x 2 ] + (-— -) , (x 2 >0) 


u 3 (x 2 } = 1 


v -,-2 


1 B l ay nR 

F-j(y) = p[- e (B-jCosB^y+BgSine-iay) " ^ * (y <0 ) > 

2v-j D-. 0 -j ay 

F 2 (y) = -pqr — e (-B 1 sin0. ) ay+B 2 cos0 1 ay) , (y>0) , 
where B^,...,B^ are obtained from the solution of the following 


_ n 9 2 "V-i 2 “Vn 

^ ^ - 2 


'I 


2E, 


) , 


b-|+b 2 = (B^-B 3 )0 2 /0i 


B 2 - -B 4 D 2 6 2 /(D 1 Bi) , 

B 2 -B 1 = (B^BgjD^/^B*) . 


(B17) 

(B18) 

(B19) 

system: 

(B20a-d) 
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Fig. 1 Geometry and notation for a reinforced cylindrical shell con- 
taining an axial crack. 


- 30 - 




Fig. 2 Normalized stress intensity factors on the outside surface of 
an "open-ended" reinforced cylindrical shell containing an 
axial through crack (Shell No. 1, R/h = 100). 
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Fig. 3 Normalized stress intensity factors on the outside surface of 

an "open-ended" reinforced cylindrical shell containing an axial 
through crack (Shell No. 3, R/h = 35.5). 
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6 



Fig. 4 Normalized stress intensity factor at the crack tip X2=0 
on the outside surface Xg=h/2 in "open-ended" reinforced 
cylindrical shells (see Table 1 for the dimensions of the 
shells and for the values of a=-y for each shell). 
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Fig. 5 Normalized stress intensity factor at the far end X£=-d 

of a through crack intersecting the interface in "open-ended" 
reinforced cylindrical shells (see insert in Fig. 4 for 
geometry and Table 1 for dimensions of the shells). 
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